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I will briefly present my work on cosmological parameters estimation. Classical methods for 
parameters estimation involve the exploration of the parameter space on a precalculated grid 
of cosmological models. Here we try to estimate the cosmological parameters by using a 
minimization method associated with the interpolation of the Ct spectrum. We first use a 
simple multidimensional linear interpolation, and show the flaws of this method. We then 
introduce a new interpolation method, based on a physical description of the location of the 
acoustic peaks in the power spectrum 



As the data available from CMB experiments improve, the estimation of the cosmological 
parameters become more and more promising. But as the power spectrum estimation increase 
in precision, we encounter new problems. The time spent calculating the theoretical power 
spectrum and the storage space are two of the main problems. Here we present a method based 
on a grid of pre calculated models with an interpolation method which compute the power 
spectrum for any given set of cosmological parameters, in order to fit the parameters using the 
most recent data available. 



1 Cosmological parameters estimation 

Our problem is a classical parameter estimation problem, but in the case of the estimation of 
cosmological parameters from the power spectrum of the CMB temperature anisotropies, the 
difficulty arises from the fact that the computation of the cosmological models is a CPU and 
storage intensive task. It takes about a minute on a standard PC to compute a single model, so 
computing the statistical quantity to be minimized (either the \ 2 or the likelihood function) in 
the process of fitting would be time consuming. One possibility is to calculate the cosmological 
models before the estimation on a discrete grid of the parameter space. Methods involving 
Monte Carlo Markov chains try to reduce the number of models to be calculated (e.g. Verde ™) 



Table 1: Cosmological parameters used for the analysis presented in this paper 
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In the case of the exploration of the parameter space, we have to compute the x 2 or the 
likelihood function for each point of the grid in order to obtain the best parameters (e.g. Bcnoit 

et 

In the case of a general minimization algorithm, we have to interpolate the power spectrum 
in order to be able to calculate the x 2 ° r the likelihood function for any given set of parameters 
ilj. This method allow us to study the degeneracies among cosmological parameters via the 
covariance matrix. 

In this analysis, we used a grid with the cosmological parameters presented in tabled This 
has a total of 640584 nodes and represents 8 Go of disk space. It was computed using CAMB 3 
with C e up to t = 3000. 

2 Cosmological parameters adjustment 

We are trying to fit the cosmological parameters, so instead of calculating the x 2 or the likelihood 
function for each point of the grid ^ we want to have a power spectrum for any given set of 
parameters fij. To realize that, we interpolate the power spectrum in the parameter space. 

We use an "HyperCube" to represent and discretize the cosmological parameter space. Each 
point of this HyperCube is defined by the cosmological parameters and the corresponding power 
spectrum. It allows us to have an irregular sampling of the parameter space. 

2.1 Simple linear interpolation 

We first use a linear multidimensional interpolation. As shown in figure ^ (left), the power 
spectrum for a given point Q{ Ci is the weighted average of the Ci{t) over p = 2 N closest 
neighbours nodes on the grid, where ./V is the parameter hyperspace dimension : 

Cintif) = X>(i)Ci(*), «>(*) = f(di,d 2 ...,d n ). (1) 
i=0 

In N dimension, the function w(i) is given by the product of ./V terms, one for each dimension 
which are either (1 — di) or (di). For example, in two dimensions w(l) = d\ X (1 — d 2 ) with 
notations according to figure ^a (left). 

As it is shown on figure|2Ja (left), this methods suffers from one major drawback. As we use 
cosmological models with different cosmological parameters, we sum power spectra with shifted 
acoustic peaks. The effect is to erase the structure of the peaks, specially at high multipoles. 

2.2 Acoustic scale 

In order to solve this problem, we use the acoustic scale defined on the figure ^ 
"http:/ /camb.info 



1 „ : sonic horizon at recombination 




Figure 1: a. left : Linear multidimensional interpolation scheme, b. right : Acoustic scale. 



With the notation of the figure, the physical size of the acoustic oscillations at decoupling 
and the angular diameter distance are given by : 
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Figure 2: CAMB computed and interpolated power spectra for £ — 800 and £ = 1600. a. left : Simple multilinear 
interpolation b. right : Improved interpolation using the acoustic scale. 

2.3 Improved interpolation 

Our new interpolation scheme is the same as on figure ^ a except that the I axis is rescaled 
according to the ratio of the corresponding acoustic scales. To compute the C(t) spectrum 
for the set of parameter {Qt} using the 2^ grid points , we define 2^ rescaling coefficients 
ai(Qt,&i) = as(Oi)/as(J7j) The interpolated C(£) spectrum can then be written as : 

2 JV 

C in t(£) = Y,w(i)Ci(lxai, w(i)= f(d 1: d 2 ...,d n ). (3) 

i=0 

Figure|21b shows the improvement of the interpolation result, compared to the simple multilinear 
interpolation. 

3 Statistical tests 

In order to quantify the improvement of the interpolation, we have computed two statistical 
quantities : 



we compute the maximum relative difference Diff max between a CAMB computed CMB 
spectrum and an interpolated one over the Ce spectrum 

we associate an error on to each point in the C(£)i n t, and we compute the Ax 2 between 
this "experiment" and the CAMB original spectrum 
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Figure |21 presents the distributions obtained for the relative difference and for the Ax 2 with 
error bars similar to the one foreseen for the Planck mission (Puget et a , with {fit} randomly 
distributed over the entire parameter hyperspace. . 
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Figure 3: Simple and improved C(£) interpolation comparison, a. left : relative difference, b. right :A% 2 for 

Planck-like error bars. 

The reduction in the mean value and the spread of the two quantities (Diff max et A^ 2 ) 
obtained by the new interpolation show clearly the improvement. 

4 Conclusion 

We have shown that our improved interpolation method provide a way to diminish the errors 
introduced in the cosmological parameters estimation process due to the parameter space sam- 
pling. Alternatively, for a given precision, it allows to use a smaller grid. It can be useful in the 
case of future experiments such as the Planck mission. Moreover, the parameter fitting is a fast 
way of estimating cosmological parameter which can be used in Monte Carlo simulation of an 
experiment. The impact of the experimental design on the cosmological parameters could then 
be estimated in this way thanks to the reduction of needed computing resources. 
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